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Abstract 

In the framework of the Hard Thermal Loop effective theory, we cal- 
culate the two-loop contributions to hard lepton pair production in a 
quark-gluon plasma. We show that the result is free of any infrared and 
coUinear singularity. We also recover the known fact that perturbation 
theory leads to integrable singularities at the location of the threshold 
for qq —> 7*. It appears that the process calculated here significantly 
enhances the rate of low mass hard dileptons. 

LAPTH-908/02, LPT-ORSAY-02/28 



1 Introduction 

Some time ago a new mechanism for the production of hard photons in a quark- 
gluon plasma was proposed [Q: the photon is produced by the annihilation of 
a quark-antiquark pair "after" one of the annihilating quark or antiquark has 
scattered in the medium (diagram on the right of figure This process is 
contained in the two-loop diagrams of the Hard Thermal Loop (HTL - 
[iH'^'Bil ) resummed theory. In the HTL approach quarks acquire a thermal 
mass in the plasma and, at lowest order in the strong coupling, the annihila- 
tion of a quark-antiquark pair into a real photon is kinematically forbidden. 
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Figure 1: The bremsstrahlung and off-shell annihilation processes that are 
expected to dominate the production of low mass hard dileptons. 

However when the annihilation is associated to scattering in the plasma one of 
the quarks becomes off-shell and the annihilation is possible. The interesting 
feature of the mechanism is that, apart from the exponential suppression factor 
common to all thermal processes, the photon production rate grows with the 
energy of the photon, leaving it as the only relevant mechanism when the photon 
energy is much larger than the temperature. In particular, it dominates over all 
other processes contained in the one- loop approximation §jl^,|ll|,^2|jl|,|l| of the 
HTL approach (qg —f q"f and qq —f gj) which only grow logarithmically with 
the photon energy. Recent studies, in the framework of the hydrodynamical 
evolution of the plasma ||l^,|l6 17 l^,|l^ , have shown that thermal production of 



a photon, including the "off-shell annihilation" process, is already significant at 
SPS energies and could dominate at RHIC and LHC energies. There is then a 
window, of photon energies of a few GeV, which could provide a signal for the 
formation of the plasma in heavy ion collisions. However, in this energy range, 
there is a very large background of pions which decay into photons. Whether 
the flux of thermal photons can be identified over the large background is under 
study. In view of this problem it may be useful to consider an alternative signal 
which follows the same dynamics as real photon production but which suffers 
from a different experimental background. The production of small mass {\/Q^) 
dilepton pairs at relatively hard energy (go) is such a process and this is what 
is considered below in the limit \/Q^ ^ qo- Considering lepton pair production 
leads to an interesting situation where two hard scales, T the temperature and 
go, and three small scales, ^/Q^, M^o ~ gT the thermal mass of fermions and 
nig ~ gT the thermal mass of the gluon enter the game. It turns out that several 
of these scales combine in a single expression Mos which acts as a cut-off to reg- 
ularize potential coUinear singularities associated with the fermion propagators. 
However, in the off-shell annihilation process, M^g can vanish and even become 
negative. Its role as a cut-off is then blurred! This occurs when the condition 

g2 > 4 Ml 

is satisfied. In this case, lepton pairs can be produced at 0(0) in the strong 
interactions by the annihilation of a qq pair, a new feature compared to the 
case of real photon production. When going to higher orders in the strong 
interactions it means that two types of diagrams will contribute, which are 
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usually referred to as real and virtual diagrams. This point, which involves the 
compensation of divergences between the two types of diagrams p0| , plp^ , |23| , p^ , 
p5[ , will be discussed in some details below as it has no equivalent in the case of 
real photon production. Even after all divergences are cancelled, the dilepton 
invariant mass spectrum still exhibits an integrable singularity at = 4 M^. 
This is in agreement with a theorem of perturbation theory which states 
that if an observable has its phase space restricted by a 6* function at a given 
order of the perturbative expansion, it may exhibit b function singularity at the 
next order. 

The collinear enhancement mechanism that makes the processes of figure |^ 
important has also been shown to play a role in multi-loop diagrams belonging to 
the class of ladder corrections and self-energy-corrections [p7|p^ , p^ , ^ , ^ 32 , 33 1 , 
and a resummation of this family of diagrams has been carried out in |32,33|. 
The effect of this resummation, also known as Landau-Pomeranchuk-Migdal 
(LPM) effect [|^j3^,^, leads to a small reduction (by about 25% for photons 
in the range interesting for phenomenology) of the photon rate. Assuming that 
the magnitude of this effect remains comparable for reasonable photon invariant 
masses (the photon mass helps to regularize the collinear singularities) , we limit 
ourselves to a 2-loop calculation. These contributions add up to the processes 
99 ~^ 91* Etnd qg — » 97* already calculated in 37 3^. 

In the following section, we set up the formalism, then we discuss the matrix 
element showing, before carrying out the phase space integration, that cancel- 
lation of some infra-red divergences can be expected between vertex and self 
energy diagrams. In Sec. 4 the phase space integration is carried out while in 
Sec. 5 the behavior of the spectrum near the singular point Q"^ — A is 
studied analytically and found to be in complete agreement with the numerical 
studies. Finally, Sec. 6 is devoted to some simple phenomenological studies. An 
appendix contains details of the phase space integration. 

In our study on dilepton pair production we rederive, as the limit when the 
dilepton mass vanishes, the results already obtained on real photon production. 
The method we follow is however simpler than previously used. It still uses the 
retarded/advanced (R/A) formalism [|9|jl|,|^,|l|,|| but it rehes on carrying out 
some integrations in the complex energy plane for which the R/A formalism is 
well suited. Furthermore, when discussing the compensation of divergences and 
calculating the remaining finite pieces, there is no need to go to n— dimensions or 
introduce regulators since the compensations can be seen within the integrand. 

Some of the results of this paper rely on a new sum-rule presented else- 
where which allows to analytically carry out otherwise difficult numerical 
integrations. 



2 Model 
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2.1 Formula for the dilepton rate 



In a quark-gluon plasma in local thermal equilibrium, the number of lepton pairs 
produced from virtual photons per unit time and per unit volume is conveniently 
expressed in terms of the imaginary part of the retarded photon polarization 
tensor, calculated in presence of the plasma Q,^: 

In this formula, the sum over the phase-space of the lepton pair in which the 
virtual photon decays has been performed, at a fixed total momentum Q = 
{ui, q) of the pair. The trace over the Lorentz indices • • • indicates that both the 
longitudinal and transverse modes of the massive photon are included. Note that 
this formula is valid at lowest order in the electromagnetic coupling constant, 
and at all orders in the strong coupling constant. In more physical terms, it 
neglects the possible reinteractions of the photon or leptons on their way out of 
the plasma, and it is valid only if the size of the system is small compared to 
the mean free path of a photon or lepton. 

Another limitation of this formula is that it holds for a plasma in thermal 
and chemical equilibrium. However, it can also be applied to a situation where 
there is only a local equilibrium, provided the typical size of the cells over which 
the system can be considered to be equilibrated is large in front of the photon 
formation timcQ |3^,Q. If this condition is not satisfied, one has to go back to 
basics, and use a non-equilibrium real-time formulation in order to compute the 
rate g7|,E8|,E9lj50|,E| R. 



2.2 Hard thermal loops 

In order to perform this calculation using thermal field theory, one must resum 
the bare propagators and vertices in order to account for the fact that medium 
effects modify the interactions and properties of the excitations of the plasma. 
In particular, medium induced masses are of utmost importance for quarks in 
processes where the photon is predominantly produced in a coUinear configura- 
tion. In thermal field theory, this is achieved through the HTL resummation, 
which resums an infinite set of one-loop thermal corrections to propagators and 
vertices. 

Since the quarks are always hard in the processes we consider throughout 
this paper, we keep only the asymptotic thermal mass from the HTL fermionic 
self-energy [ |5^ , which leads to the following retarded and advanced propagators 

S^. (P) ^ Ps,,. (P) , (2) 

^ This is the imphcit assumption in the hydrodynamical approach used in order to describe 
photon emission during the collision of two heavy ions. 

^Although the assumption of an infinitely fast switching of the coupling constant that has 
been made in thosRjcalculations may lead to unphysical effects like a power-like tail in the 
emission spectrum p2| . 
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with 



P EE (po, Vp' + MIp) 

where = g^C^T'^ /A is the square of the thermal mass of a hard quark 
{Cp = {N^ — l)/2Nc is the Casimir in the fundamental representation of the 
gauge group SU{Nc)). Ordinarily, one does not need effective jqq vertices if the 
quarks are hard, since the HTL corrections are suppressed. But this is correct 
only if the terms we are calculating are at leading order. Photon production 
however involves an important cancellation between the components IIoo and 
Hzz of the photon polarization tensor, and the remaining terms in 11^'' are 
non leading, and are affected by the HTL corrections to the vertex"^. In the 
following, we do not need the explicit form of the vertex, but only need to know 
how it is related to the self-energy correction by means of Ward identities. In 
other words, it is important to make sure that we do for the effective vertex an 
approximation which is consistent with the approximation of Eq. (|^) made for 
the quark propagator. In particular, we want to preserve the relations satisfied 
by the full HTL corrections (for a vertex in which a quark enters with momentum 
R and goes out with momentum P): 

p^sr^ - m = H-H, 

R^ST^^ = m = f'-f, (4) 

where ST is the HTL correction to the vertex. In the following, it will be useful 
to write 

^ Y^'•'-f^ . (5) 

In addition, we need the effective retarded and advanced propagators of a 
gluon, which can be decomposed into its transverse and longitudinal compo- 
nents: 

- D>^';jL) EE P^%L)^l jL) + P^'^iL)Al jL) , (6) 

where P^"^ (L) are respectively the transverse and longitudinal projector^ for 
a gauge boson of momentum L [54|| : 



•^This vertex correction has never been considered in the existing calculations of photon 
production. As we shall see later, it slightly modifies the value of 11^^'. However, one can 
check that it modif ies only the IIzz component, and therefore forgetting it had no impact on 
the result of |32y33 who calculated only the transverse components (since IIoo = ^zz for real 
photons). For virtual photons, one has to include the temporal and longitudinal components 
as well, and one cannot avoid using this vertex correction. 

*They satisfy the following useful identities: 

Pj, Pj, -Pp , P^ Pj^ P^ , P^ P^ P^ Prp , 

L^pi:-; (L) = . (7) 
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(8) 



where 7 



gfif „ U^U", = j^^Li^, and where U is the mean velocity of 



the plasma in the frame under consideration. The propagators of the transverse 
and longitudinal modes are then given by 



R,A ^ ' 

with the following HTL self-energies 



X- 1 



(9) 



(10) 

where we denote x = Iq/1 and where rrig = g'^T'^[Nc + Np/2]/9 is the gluon 
thermal mass in an SU (Nc) gauge theory with flavors. Note that an implicit 
continuation Iq —^ Iq ± ie is understood in this self-energy in order to obtain its 
retarded and advanced components. 



3 Matrix element 
3.1 Diagrams 

Let us now move to the two-loop diagrams calculated in this paper. They are 
the same as those already considered in [||J|,|l| for real photons and in ||5^ for 
soft dileptons produced at rest in the plasma frame. One of them is a vertex 
correction to the one-loop photon self-energy, and the other is a self-energy 
insertion on one of the quarks lines (see figure ||). 

There is however one major difference compared to the case of real photons: 
as long as the invariant mass of the lepton pair is < iM^, the cuts denoted 
by (a) and (c) in figure ^ do not contribute. Indeed, they correspond to the 
interference between the tree-level process qq — > 7* and a virtual correction to 
this process, but this direct production mechanism has a threshold at = 
4M^. Above the threshold, one must include these cuts, and in fact the various 
cuts partly compensate each other through the KLN mechanism. 



3.2 Vertex contribution 

It is convenient to use the retarded-advanced formalism [^9|jl^,^,El j4|] in order 
to compute the retarded photon polarization tensor needed in Eq. (n ) . Through- 
out this paper, we assume that the gluon of momentum L (see diagrams and 
notations in figure ^ is soft, so that the corresponding Bose-Einstein factor 
ng{lo) is large, in particular large compared to fermionic statistical factors. 
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Figure 2: The two- loop diagrams contributing to the production of hard dilep- 
tons. A third diagram with a self-energy insertion on the upper quark line is 
not shown here. 



When performing the integrals, one can check the validity of this hypothesis. 
Therefore, we keep only those terms containing this enhanced Bosc-Einstein 
factor, and we approximate 71,^(^0) ~ T/Iq. 

Using this approximation, the contribution of the first diagram (which we 
designate by "vertex correction" in the following) of figure ^ to the photon 
polarization tensor is given by: 

1 2 „ f d^P 



Imn4(g) « T;^'g'N,C, / — — [n,(ro) - n,(po)] 

vertex Z J [ZTT ) 



X 



a=T,L 

X [,s, iP)s„ (P + L)- {P)s, {P + L)] 
x[s,{R)s,{R + L) - s,{R)s,{R + L)] , (11) 

where pa {L) is the spectral function for the mode a of the effective gluon prop- 
agator: 

p„(L)^A^(L)-A^(i), (12) 
and where Trace™"^* denotes the trace of the Dirac matrices for this diagram: 



Trace™;* ^ Tr + <5r^)/?(7p + Srp){fi + + ST^){p + + . 

(13) 

This trace contains four hard momenta, and is a priori of order or larger. It 
is however important to keep also subleading terms of order T^M^ since the 
cancellation between IIoo and kills all the leading terms when the vertex and 
self-energy diagrams are combined. However, it is safe to drop all the smaller 
corrections, i.e. terms of order or smaller. Expanding this trace up to 
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terms with two powers of a soft quantity, we find^: 
-4 



Trace 



per 



2(i? +(i? + L)^)PpQ,-2(P +(P + L)^)i?pQ, 
-8(P" + T"^Pp)(Ro. + TorR^,)PpRa 



r"^Qp){Qa+Ta^Q^)PpRa 



(14) 



In this expression, one can even drop the terms with two powers of the vertex 
correction F'^" , since they would be smaller than the terms we are interested in 
here. Using the pseudo Ward identities of Eqs. (^), it is possible to rewrite the 
three terms in PpR^ as 

16 [Q^ + 2{P ■ R- P ■ R)] PpR, , (15) 

up to terms that contain more powers of soft factors. 

3.3 Self-energy contribution 

Similarly, the contribution of the diagram with a self-energy insertion on the 
quark of momentum R is given by 



Imn^'KQ) 



1 



self 



(PL T 



d*P 
(2^ 



[npiro] - n^(po)] 



(27r)4 ^0 



^ p„(i)Pr (i)Trace^pf 



=T,L 



x[[s,{R)]h,iR + L)-[sAR)?sAR + L)] , (16) 



where Trace^^^ is the corresponding Dirac trace: 

TVace^f = Tr [^(7^ + SV^ jfiijp + ^r,)(;^)(7, + STM^/p + STp)] . (17) 
Calculating this trace down to the order T^M^, we obtain 



Trace"; 



AlfRpQa + 4(P^ + lf)RpRa 

-8(P" + r"'5p0)(P„ + T^^R^)RpR, 

-8(Q« + r"^Q^)(Q„ + r„^g^)PpP, 



(18) 



We can recognize the same set of terms as the one evaluated in Eq. ( |l5| ) , with 
PpRa replaced by RpRa- Of course, there is also another such diagram, not 

^Anticipating the fact that contracted with the projectors P^*^ {L) gives a vanishing 
result, we have dropped the terms proportional to this quantity. 
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represented explicitly in figure g, with the self-energy insertion on the quark 
of momentum P. In fact, the contribution of this diagram is equal to that of 
the first self-energy insertion. Alternatively, the corresponding matrix element 
can be obtained from the previous one by performing the following change of 
variables: 

-P , 

L^-L. (19) 



This second method, compared to the one where we just multiply Eq. ( |18| ) by a 
factor 2, has the advantage of making the matrix element symmetric under the 
exchange of P and R and is preferred in order to see some cancellations. 

3.4 Partial cancellation between vertex and self 

Without performing any explicit calculation, one can first prove that the terms 
in RpQa and PpQa all cancel when one combines the contribution of the vertex 
correction diagram, and of the two self-energy insertion diagrams. These terms 
can be rewritten as^ 



/-r:2 



2{R +{R + LY)PpQa {P)s^ {P + L)s^ {R)s^^ [R + L) 



-2(P + (P + LY)RpQ, s^{P)s^{P + L)s, {R)s, (R + L) 
+m^RpQ, s^{P)[s^{R)Ys^{R + L) 

-iP^PpQa s^iR)[s^{P)fs^iP + L)\ +C.C. . (20) 

2 

Using P a(P) = i, one can check that the square bracket in this equation is 
proportional to: 

'aPpQ, (P)s, (P + L) (,s, (R + L)- (P)^ 

~4RpQ^ sAR)sAR + L)(s„{P + L) - s,{P))] . (21) 

Assuming that the gluon momentum L is small compared to the quark momenta 
P and P, and performing the change of variables 

P + L^-P, 

i ^ P , (22) 

^The complex poles of the statistical facto rs d o not crmtribute since the difference between 
retarded and advanced propagators in Eqs (hll) and (llm vanishes when evaluated at such 
poles. 
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on some of the terms (dropping terms in Lp that may appear), one can observe 
that the two lines cancel each other. This cancellation is in fact a particular case 
of a cancellation between ladder diagrams and self-energy insertions discussed 
in ||2^,||,|6|j5^,|5|,||,|2|Q. For the remaining terms, we have: 



Imn,;i(Q)«-2eVA^cC 
T 



K(ro) 



(2^)4 ^0 



a=T,L 



2L^{RpRa+PpPa) s^{P)s^{P + L)s^{R)s^{R + L) 
-i{Q^ + 2{P-R-P- R))RpPa {P)s^ [P + L)s^ {R)s^ [R + L) 
-A{Q'' + 2{P-R-P- R))RpR, s^{P)[sa {R)?s^{R + L) 
-4(Q2 + 2{P-R-P- R))PpP, ■^A{R)[sn{P)?s^{P + L]^ 

(23) 



Let us add that, since we are looking at photons whose virtuality is not negli- 
gible, we have kept the terms in when they have the coUinearly enhanced 
quark propagators. 



4 Integration over the quark momentum 

4.1 Pole structure of the integrand and collinear enhance- 
ment 

We are now in a position to perform the integration over the quark momentum. 
We perform this integration in the complex plane using the theorem of residues, 
and we integrate first over the component of the quark momentum which is 
parallel to the photon momentum^ i.e. the z-component. 

For each term, we have propagators with a retarded prescription multiplied 
by propagators with an advanced prescription. It is convenient to perform in 
the complex plane the integration over the longitudinal components of the quark 
momenta (the reference direction is the photon momentum) , by closing the con- 
tour around the upper half-plane. In this approach, the collinear enhancement 
is due to the fact that in a product like s^(P)s^(i?), there are pairs of poles that 
are separated by a very small interval and are on opposite sides of the real axis. 
As a consequence, the real axis is "pinched" by such a pair of poles, leading to 
a large contribution of the order of the small separation between the poles. 

Keeping only the contribution of the pinching pole, Pz = po+{p\+M'^ff) /2po, 

^We chose the photon momentum Q to be along the z axis. Therefore, by definition, we 
have = 0. 
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the basic result we need is the following (for R = P + Q) 

/■+°° dv 11 
J ^ i^'n{P)sAR)F{p.) « ^^^^^72 71^F{p,) , (24) 

where, foUowing |§,|l|, we denote 



Af2 _ jBO 

cfi go 



M^s = Ml + ^poro . (25) 
% 

This resuh is valid if the function F is not participating in the colhnear enhance- 
ment (i.e. is well behaved in the vicinity of the poles of the quark propagator). 

4.2 Contribution of the vertex diagram 



For the calculation of the first two terms in Eq. ( p3| ) (i.e. the terms coming 
from the vertex correction diagram), we use this method in order to perform 
the integral over the z-component of the quark momentum in each loop, i.e. 
respectively pz and pz + Iz considered as independent variables. These terms 
become: 

imn^^:(g)| ^ « 

H^\^ / \ vertex 

eV^cC^ [ dpo(fp^ fdlodH^T 



J2 pa(i)pr(i) 

[L^iRpRa + PpPa) - 2(0^ + 2{P-R-P- R))RpP, 



X 

Q=T,L 



Pz=PO 



1 1 



" Pi + - (p^ M^^ - ^I^e + ' ^^^^ 

Next, we need to perform the contractions with the projectors. For any pair of 
4- vectors A^, taken in the set {Pp,, -R^}, we have in the coUinear limit: 

A^B.P^'^iL) - (i ■ a){i • 6) - (a • 6) « -A^B^P^^iL) , (27) 

where I = I /I. Noticing that all the transverse components of p and r are much 
smaller than the longitudinal components due to the coUinear enhancement, we 
have also ^ ^ 

i-P-Po-y^ i-rfuro-^, p ■ r ^ pqTo , (28) 



so that 



I 

,12 
2 

,1-2 



PpP.P^^iL) « -PpP,PP^%L) ^pli^-^-lY 
RpR^P^iL) « -RpR,PP%L) « - l) , 

PpRaP^" {L) « -PpR^PP^ (L) « poro (| - l) ■ (29) 
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It is convenient to use the variable x = Iq/I instead of Iq itself. However, since 
Iz = Iq and P = Ig + l']^, the Jacobian of the transformation differs from 1. A 
simple calculation gives: 



(30) 



Therefore, this change of variable enables one to absorb all the factors 1~ Iq/P 
coming from the contraction with the projectors*. 

Finally, taking m = and v = l'^ as the integration variables, we can 
rewrite after carrying out the angular integrations: 



imn,j:(g)| 



-f OD 



1 +00 

^ / dv 



X 



{v + 2Ml)ipl+rl) + 2Q^Poro 











2Imn^(a;) 



2Imn^ {x) 



(w + Ren^(a;))2 + (Imn^(a;))2 {v + Ren^{x))^ + (Imn^(x))^ 



du- 



sign('u 



(u + V + Af2jj)2 - A.UV 



1/2 



(31) 



where we denote 



M 



2 — n/r2 

off 



Poro 
f 

90 



(32) 



This ie prescription coming from the quark propagators will turn out to be 
important when M^^ < 0. 

4.3 Contribution of the self-energy diagrams 

A similar strategy can be pursued for the contributions of the self-energy dia- 
grams, i.e. for the third and fourth terms in Eq. (p3|). Using the fact that these 
two terms are equal, and performing the integration over pz, we first obtain 



self 
2„2 



Ae^g^NcCj, f dpod^p 



(2^) 



d'^L T 



X [p^ [L) - p^{L)] {Q' + 2{P ■ R- P ■ R))r\l - cos' 6ri) 



* Since L is space-like (as required for scattering processes like those of figure ih, the inte- 
gration range for the variable x is [— f , +1], which can be replaced by [0, 1] using the parity of 
the integrand in the variable x. 
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(33) 



where 9ri is the angle between the vectors r and I. The integration over this 
angle is straightforward, since we have: 



+1 



d cos 9- 



rl 



TqIo — ir^e — rl cos dri 



1 

Vl 
1 

Vl 



1 - 



(»^o^o) 



21 



{rlf 



In 



^^0^0 - irpd + rl 
roh - iroe - rl 



sign(ro) 


In 




+ in 


lo-l 









rplp 
{rlf 

rplo 

'[riy 



(34) 



Noticing that for soft the factor p{L)T/Iq is an even function of lo, we can 
drop any term in the above result which is odd in Iq. Therefore, only the imag- 
inary part contributes. Using this result, and introducing again the variables 
X = lo/l, u = p\ and v = —I? = P(l — a;^), we can write: 



327r4 



1 +00 





X 



2g%ro + 2M^(p2+rg) 
2Imn^(a;) 



2Imn^ (a;) 



{v + KeIV^{x)Y + (Imn^(a;))2 (w + Rcnja;))^ 
X j du- ^ — + c.c. , 





(Imn,(x))2 



(35) 



4.4 Cancellation between vertex and self 

Combining Eqs. (|3^) and (|5|), one obtains 



inin,i:(0) 



+ 00 



oo 4 -2 1 dpa[n^{ro)-n^{pQ)\ — dv 

OZTT Qo J J X J 

-oo 



dx 



+ 00 



2Imn^(a;) 



2Imn^ (a;) 



{v + Ren^{x)f + (Imn^(a;))2 {v + RcIl^{x)f + (IrnHja;))^ 

+ 00 



du 



u + M^^ 



sign(M + V + M^g) 

1/2 



{u + v + Af2jj)2 - iuv 
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sign(u + V + M^g) 



1/2 



ctt 



+C.C. 



(36) 



Therefore, we see that there is a canceUation for the term in between the 
contribution of the vertex and the contribution of the self-energy corrections, 
in the hmit where w — > 0, i.e. when the momentum transferred by the gluon is 
small. This is an extension to virtual photons of a well known cancellation in 
the case of real photons |27||2^ , |5^ , p7 58 2^, |3^ , p3[ . In particular, it ensures that 
the dilepton rate is not sensitive to the gluon magnetic mass if the magnetic 
mass is small enough (say of order g^T). 



4.5 Integration over p\ 

At this stage, it is possible to perform analytically the integration over the 
transverse momentum of the quark u — p\. One has a priori to distinguish 
two cases: M^g > and M^^ < 0. The simplest case is M^g > 0, for which 
the denominators never vanish and for which the ie prescription is irrelevant 
(therefore, adding the complex conjugate just amounts to multiply by a factor 
2). In this case, we have: 



+ 00 



du 



sign(w + V + M^g) 



8tanh"^ 





+ C30 



du 



(u + v + Afcff )^ - 4:uv 
2 



1/2 



c.c. = 



-0+ M„^ 



0+ ^v{v + 4Af^g) 

(37) 



A similar calculation can be carried out when M^g < 0. Note that in this 
case the result of the integral over u is a complex number, but because we 
need to add its complex conjugate, only the real part is important to us. In 
addition, the calculation shows that we have to distinguish according to whether 
V + 4Af^jj > OT v + 4Mgg < 0. The results are the following (see the appendix 
^ for details): 



• If 

4-00 



M^ff < and v + AM^g > 
du sign(w + v + M^g) 



Stanh 



u 



• If 

4-00 



1/2 



{u + v + M^f^y - iuv 
M^ff < and v + AM^g < 

du sign(M + V + M^ff) 



-0+ ^v{v + ur^^) 



Stan 



M, 



off 



{u + v + M^ff Y - 4:uv 



1/2 



-0+ ^-v{v + 4M^g) 
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(38) 



One notices that those results could have been obtained easily by taking the 
real part of the analytic continuation of the M^g > result. Indeed, if M^g < 
and V + 4Mgfj > 0, then \Jv/{v + 4Mgfj) > 1 and this analytic continuation 
just amounts to replace the argument of the tanh^^ in Eq. ( |37| ) by its inverse. 
If on the contrary v + 4M^fj < then the square root in Eq. is a purely 
imaginary number, and by using iXtmi\i-^{iX) = -Xtan-i(X) one obtains 
easily the correct answer. 

Denoting v = M^f^w and introducing the following set of functions^ 

+00 



\/w/(w + 4)tanh + 4) 

F2 



=^^'^/T^"^"--(^)/'^^(M> + Ren,,(a:))2 + (Imn,,(x))2 ' 

00 ■ ' 

1 +00 



1 +00 _^ 

dx ^ ^ ^ ^ f dw y/ w / {w + 4:)ta.nh y/w/ {w + 4) — w/4: 

^(M2gu; + Ren^_Jx))2 + (Imn^_Ja;))2 ' 

(39) 







one can write a very compact expression for the 2-loop photon polarization 
tensor: 

+00 

Imn4(g)« Vr^- / dPoKiro) - n^Po)] 



Ml 



y-{{pi + rt){J^ -J.) + 2 ""'^" 'Zr" "' (K. - K.) 

off 



(40) 



Note that the above defined ^ match those already defined in ||^ . Having this 
in mind, one can recover the limit of real photons. Indeed, in the limit where 
^ the quantities and ^ become independent of the quark energy 
Po, and can be factored out of the integral (they depend only on M^jm,^. 
Therefore, we obtain: 



lmImn,j;(Q) « \_l - -[J^-J^+2K,-2K,) 



X j dpQ[n^{ro) - n^(po)](Po +''0) 

— 00 



-(J^ - +2K^ -2K^ 



2^3 

qoT + n^ — 

10 , 



(41) 



'Or their analytic continuations if A/|?fj is negative. 
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Note that the integration over po is performed exactly her^^ In the papers 
the numerical evaluation of J^^ was overestimated by a factor of 4, as 
first pointed out in and independently in [Q. We also see that taking 

correctly into account the HTL correction to the vertices, which was ignored in 
[Hl^, brings another term of the same order, namely 2K^ — 2K^. Numerically, 
this term is a 30% negative correction for Nc = 2> and = 2 or 3. Although 
the references |^,^ did not consider this HTL vertex correction, they obtained 
the correct answeiP] for real photonsQ because they chose to calculate only the 
transverse part of the polarization tensor (XIn +n22)- Indeed, it is easy to check 
that the HTL correction to the vertex modifies only the component lizz ■ 

One can also add that Eq. (|4n) combines in a single integral the contributions 
of all the processes in figure ^ The bremsstrahlung is obtained for po > 
(bremsstrahlung of a quark) and po < ~<lo (bremsstrahlung of an antiquark), 
while for —go < po < one gets the off-shell annihilation. 



4.6 Exact integration over x 

The advantage of the expression given in Eq. (^) is that the integrals ^ and 
Krj. ^ it contains are functions of Af^g /mg which can be studied rather simply 
analytically. Indeed, we show in a separate paper |Q that the result of the 
integration over the variable x is in fact extremely simple thanks to the use 
of sum rules. We show that for a general enough self-energy (see |Q for the 
conditions under which it is true), one has: 



dx 



2ImH(2:) 



X (z-f ReH(a;))2 + (ImH(a;))2 



1 



1 



^ + ReH(cx)) z + ReH(0) 



(43) 

From there, it is possible to give integral expressions for the functions ^ 
and K^^ that are much simpler than Eq. (^ ) and are very well suited for 
the derivation of various asymptotic expressions. In this paper, we just quote 
the results we need without proof, and refer the reader to |Q for a thorough 
justification. For instance, we obtain: 



du- 



tanh ^ {u) 
8 V3m2 



-'Indeed, one can prove: 



dx 



1 



3TOg 



du 



tanh ^ (u) 

(^^-1)^2, 



1 



-I- 1 e^+y + 1 



.(44) 



(42) 



^^But we do not understand how they could find an agreement with our incomplete result 
of lyl, as stated on page 8 of |B3]. 

^^heir result is incomplete lor virtual photons, because they did not include the longitu- 
dinal mode of the massive photon. 
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For real photons (Q^ = 0, Af^fj = and for Nc = 3 cofors, we have: 

4A//2„ 8 



3m2 6 + iV^ 



(45) 



Therefore, there is an accidental simplification for = 2 flavors, as this ratio 
is then equal to 1. In that case, we have 



K, -K. 



TT ln(2) 



= -(l-21n(2)) , 



and therefore the 4-dimensional integral of Eq. (p6| ) can be done exactly: 



Jim luiU^^JQ) 



2„2 



Nc=3 
N„=2 



2e^ 
' 37r3 



'11 



(46) 



(47) 



Note that this expression is valid for hypothetical quarks of electrical charge 
e. If the two quark species under consideration are u and d quarks, it must be 
multiplied by a factor (1/3)2 + (2/3)^ = 5/9. 

5 Behavior near the tree-level threshold 



5.1 Analytic expressions near Q'^ 



4M2 



A close look at Eq. ( ^0[ ) seem to indicate that there could be problems due to 
the denominator M^f^ for the second term since this effective mass parameter 
can vanish if > 4M^. This is not a problem as long as the zeros of M^g (in 
the variable pq) are simple zeros, since this term should be understood with a 
principal value prescription. This is the generic case, since the zeros are simple 
for any not equal to 4M^. However, there is a problem for a photon invariant 
mass — 4:M^, i.e. at the threshold for the tree level process qq^ ^* ■ Indeed, 
for this value of Q^, the quantity Af^g has a double pole at po — ~qo/'^- Being 
a double pole, it cannot be dealt with a principal part prescription and makes 
the result infinite. 

In order to make this statement more precise, we derive here the analytic 
behavior of Imn^^((5) for in the vicinity of 4Af^. This calculation is made 
simple by the fact that we want to extract only the diverging pieces near = 
4:M^. We can therefore drop the term in since it is finite, and replace 
Po^oQ^ by — gpAfj^ (since « 4Af^ and po « —qo/2 is the location of the 
singularity), which gives 



imn 



threshold 



27r4 



+ 00 



■TMi 



dpo[np{ro) - n^(po)] 



Af2„ 

ctt 



(48) 
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In order to simplify further the calculation, let us limit ourselves to a large 
photon energy qo » T. In this case, one can check that the statistical weight 
'^F (Po) ~ (Po + Qo) can be approximated by a function which is 1 in the range 
[—90,0] and elsewhere. Doing so is accurate up to terms suppressed by at 
least one power of T/qq ^ 1. Due to a symmetry of the integrand, we can in 
fact integrate only over [—go/2,0] and multiply the result by a factor 2. 



5.2 Case < 4M^ 

At this point, we have to distinguish two cases, depending on whether < 
4:M^ or > 4M^. Let us start with the case < 4M^. In this case, it 
is convenient to replace the integration variable po by a variable z defined as 
follows 

qo ^ qgMo. 
2 ~ ^ 



Po 



'1 - 



4M2 



(49) 



The range po E [-qo/2,0] is mapped onto z e [0, a/Q2/(4M^ - Q^)] in this 
transformation^ and the expression of M^g becomes: 



M^g = -{AMl~Q'){z' + l) 



(50) 



This leads to the following expression for the photon polarization tensor near 
the threshold: 



Imn4(Q) 



qoTMl 



+ 00 



:=(4M^)- 



VQ2(4M^-Q2 



dz ■ 



z2 + 1 



(51) 



At this point, we need only to know the behavior at small positive M^g of the 



functions 



In 



we prove the following results: 



M^ff<^mg 8 



In 



Clt 



7rln(3) 



(52) 



where the terms neglected vanish when M^g — > and therefore do not contribute 
to the singular behavior near the threshold. The longitudinal contribution to 
this self-energy is therefore given by 



d(4M2 



87r2 



VQ2(4Ae-Q2) 



(53) 



^■'Note that the upper bound for z is large near the threshold, and can be replaced by +00 
in this calculation. 
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while the transverse contribution is^ 



In 



(55) 

Note that these terms arise only from the cuts (b) and (6') of figure ^ (the cuts 
(a) and (c) are zero in this domain of Q^). 



5.3 Case > AM^ 

If the invariant mass squared is above the threshold, we need to make a different 
change of variables, because M"^^ has zeros in the integration range. Now, the 
appropriate change of variables is 



, ^0 



AMI 



l2, 



which leads to 



(56) 



(57) 



Again, when we are very close to the threshold, the range [—go/2, 0] is mapped 
on z g [0, +oo[, which enables to write the following expression for the photon 
polarization tensor above the threshold: 



Imn4(Q) 



-f2 



dz ■ 



Z2- 1 



(58) 



The pole at z = 1 should be handled with a principal value prescription. We 
can make use again of Eqs. ( p2| ) (with the difference that the absolute value 
of Mgjj should enter in the logarithm for K^), and obtain for the longitudinal 
contribution: 



imn,j;(Q) 



d(4M?,)^ 



4-^-3 



VQ2(g2 -4Ae) Jo z 



dz 



1 







(59) 

However, the integral over z is vanishing due to the principal part prescrip- 
tion. In other words, the longitudinal contribution has no singular part above 
the threshold. For the same reason, the non vanishing singular piece in the 
transverse contribution is 



imn,{:(g) 



qoTM^ 



47r3 



VQ2(Q2_4M^)7o 



, lull 
dz — : — 



^^For this calculation, we need the following definite integral: 

f 

Jo 



ln(z^ + 1) 
dz \ ' = 7rln(2) 



+ 1 



(60) 
(54) 
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In order to compute the final integral, it is convenient to split it in two terms 
z £ [0, 1] and z G [1, +oo[, and to perform the change of variables z — > l/z in the 
second term. Upon applying this transformation, the integral over z becomes 

^ ln|l-z2| .1 Ij^(^) ^2 

Therefore, the singular part of the transverse contribution above the threshold 

The slightly less singular behavior obtained for > 4:M^ is due to a par- 
tial cancellation between the cuts (a), (c) and (6), (b'). One can also note that 
the singularities exhibited in this section are integrable: the energy spectrum, 
integrated over the photon mass, is therefore always finite. 



5.4 Numerical illustration 

The above formulas for the singular behavior of the photon rate near = 
4:M^ can be checked by a full numerical evaluation of Eq. (j^). The numerical 
evaluation is done for Nc — 3 colors and Nf = 2 light flavors, with a strong 
coupling constant g — 2 (i.e. « 0.32) and a photon energy qo/T — 50. 
In figure ||, one can sec clearly that there is no divergence at — 4:M^ in 



1000 




0.01 I ■ ■ — ^ ' — ^ ■ ' — ^ ■ — i ' ' — ^ 

le-05 0.0001 0.001 0.01 0.1 1 

Figure 3: Behavior of the transverse and longitudinal 2-loop self-energy just 
above the tree-level threshold. The behavior of the transverse contribution is 
compared to the analytic formula Eq. ( |62| ) for the singular terms. Dimensionful 
quantities are in units of the temperature T. 
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the longitudinal contribution, and that there is a divergence in the transverse 
contribution, correctly reproduced by Eq. (p^). 

Figure ^ presents the same results below the threshold. Here, both the 
transverse and longitudinal contributions are singular, in very good agreement 
with the analytic prediction of Eqs. ( p5| ) and (|5^). 

10000 r ' ^ • • ' ^ • 




le-05 0.0001 0.001 0.01 0.1 1 

Figure 4: Behavior of the transverse and longitudinal 2-loop self-energy just 
below the tree-level threshold. The numerical results are compared to the ana- 
lytic formulas Eqs. (|5|) and (|5^) for the singular terms. Dimensionful quantities 
are in units of the temperature T. Note that we display — Imllj^ and not Imll^. 

We observe that the 2-loop photon polarization tensor considered in this 
paper is not defined if = 4M^, and more generally that the perturbative 
expansion leads to inconsistent results in the vicinity of the tree level threshold. 
This could have been expected on very general grounds |p6| . Indeed, it is well 
known in perturbation theory that if a contribution at order n has a phase- 
space constraint (like the 0{Q^ — 4M^) for qq 7* at tree level), then higher 
order corrections to this contribution exhibit a singularity at the point where the 
constraint starts. A correct assessment of the behavior of this process near the 
phase-space boundary requires usually the resummation of an infinite number 
of terms. 

In our case, among the next order corrections are a correction SM^ to the 
thermal mass of the hard quark, as well as a width F for the quark. Particularly 
important is the width which must be resummed whenever M^g is small [ p9| , 
(and this threshold problem is due to the possibility that M^g vanishes). 
However, gauge invariance dictates that ladder corrections to the qq^ vertex be 
also resummed. Therefore, one can anticipate that a complete treatment near 
the threshold involves a simultaneous resummation of a width on the quark 
propagator, and of the ladder corrections to the vertex where the photon is 
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attached. Since virtual photons have a physical longitudinal mode, doing this 
requires to extend the work of (3^:^ in order to resum also the LPM corrections 
for the longitudinal mode p^ . 

5.5 Positiveness of the rate 

Away from the threshold, the 2-loop photon polarization tensor if finite, but 
may have the wrong sign in order to give a positive dilepton production rate, 
as illustrated in figure pi We see that this self-energy would lead to a negative 





-0.5 



-1 I ' — ■ — ' ■ — ■ — ■ ■ — ' — 

0.1 1 10 100 

Figure 5: Sum of the transverse and longitudinal contributions to the 2-loop 
photon polarization tensor. The photon energy is set to qo/T = 50. We see 
clearly that there is a region of where —lmIl^{Q) is negative. Dimensionful 
quantities are in units of T. 

contribution to the photon rate immediately above = 4M^. This by itself 
is not enough to indicate a violation of unitarity since above the threshold the 
2-loop contributions are nothing but higher order contributions to the tree-level 
process. They may be negative, but are suppressed by a power of , so that 
the total rate should remain positive. 

Let us illustrate more graphically this point. The 1-loop diagram contains 
only the direct production of a virtual photon by the annihilation of a quark 
and an antiquark 




However, when we consider the following two-loop diagrams, we have not only 
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new processes like bremsstrahlung and ofF-shcll annihilation 



+ 



but also the interference between the Born level amplitude and loop corrections 
to them 



+ 




c.c. 



Naturally, these interference terms can be negative, and are responsible for the 
fact that the 2-loop contribution is sometimes negative. Therefore, even if we 
add up the tree level contribution 



imn,;i(g) 



Born 



An 



(63) 



it may happen that we get a negative rate in the vicinity of the tree-level thresh- 
old, due to the breakdown of perturbation theory at this particular point. This 
is in fact what happens, as shown in figure ^. We show on this plot that the sum 
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Figure 6: Born term, 2^2 processes, and 2-loop contributions. Also shown is 
the total yield (up to 2-loop) with a finite mass resolution. The photon energy 
is set to qO/T = 50. Dimensionful quantities are in units of T. 

of the Born term and the 2-loop terms has the wrong sign just above ~ 4M^. 
This indicates that the 2-loop result should not be trusted in this region. The 
solid curve shows that by dividing the range in finite size bins (as would be 
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the case in any realistic experimental situation, where there is only a limited 
mass resolution), one may average out positive and negative contributions, and 
get a rate that is always positive. This observation should not be considered as 
the definitive cure of this problem, which should consist in a resummation to all 
orders. Also shown on this plot for the sake of completeness is the contribution 
of the processes qq g^f* and qg qj* , calculated in Q . These processes 
appear in the 1-loop HTL photon polarization tensor, and contribute also at 
the order e^g^. The corresponding contribution to Imll^^ is given by: 

(64) 

Note that this expression should not be extrapolated to very small photon 
masses. At small Q^, the in the logarithm is eventually replaced by the 
quark thermal mass [ p^pT|. In order to take simply this effect into account, we 
limit the growth of Eq. (|64|) at small by replacing it by the ~ result 
of whenever Eq. ( p^ would give a larger result. It would be interesting 

to recalculate these 2^2 processes for photon masses comparable to thermal 
masses in order to obtain a more correct matching between Eq. ( |64| ) and the 
Q2 ^ hmit. 



Imn ''(O) eg ^cC^ 2 



6 Phenomenology 

In this section, we present some results in a less academic situation. We choose 
parameters as they may appear for heavy ion collisions at LHC. The temperature 
is set to 1 GeVp|. We take a coupling constant g = 2, i.e. ~ 0.32. The 
number of colors is set to Nc — 3 and we take two flavors (u and d, with 
respective electric charges 2/3 and —1/3). 

6.1 Mass spectrum 

One can plot first the mass spectrum, at a fixed photon energy. In figure 0, the 
photon energy is set to go = 5 GeV, and we plot the yields for masses between 
200 MeV and 3 GeV. It appears that at this value of the energy, the 2-loop 
processes we compute in this paper is comparable or even slightly larger than 
the 1-loop HTL contribution, especially in the vicinity of the threshold. At 
the same energy of 5 GeV, this process is slightly more important for a lower 
temperature, as illustrated in figure ^ 

The other important remark is that at such values of the coupling constant, 
the thermal masses are not small, and the threshold for the Born process is 
located at rather high masses. Therefore, most of the spectrum is in fact dom- 
inated by formally higher order terms. 

^^This may be on the high side, but most of the photons and light dileptons are produced 
in the early stages of the plasma evolution, when the temperature is still very high. 
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Figure 7: The dilepton rate in physical units, for Np = 2, T = 1 GeV, g 
and a fixed energy go = 5 GeV, as a function of the mass of the pair. 
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Figure 8: The dilepton rate in physical units, for = 2, g = 2 and a fixed 
energy qq = 5 GeV, as a function of the mass of the pair, for two values of 
the temperature T = 1 GeV and T = 0.5 GeV. The vertical lines indicate the 
contribution of the terms calculated in the present paper to the rate. 
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Figure 9: The dilepton rate in physical units, for ~ 2, T — 1 GeV, g = 2 
and at two fixed pair masses {Q = 400 MeV and Q = 1.5 GeV), as a function 
of the energy of the pair. AU the contributions of figure |^ are added up. 

6.2 Energy spectrum 

In figure ^, we set the photon mass to some fixed vakie {Q = 400 MeV and 
Q = 1.5 GeV), and we plot the dilepton yield as a function of the energy, for 
energies between qo = and 10 GeV. Note that in the region where go ~ Q) 

the approximations made in this paper as well as in |^ for the 2^2 processes 
are a priori not valid. We observe a very fast drop of the yield with energy, 
following the usual exponential law in cxp^—qo/T). 
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7 Conclusions 

In this paper, we have calculated the contribution of bremsstrahlung and off- 
shell annihilation to the production of a high energy low mass dilepton by a 
quark giuon plasma. As long as the photon mass remains small in front of its 
energy, this process is coUinearly enhanced in the same way as for the production 
of a real photon. Cancellations between real and virtual cuts ensure that all 
infrared or coUinear divergences cancel. 

Interestingly, our result display a general feature of perturbative expansions: 
such an expansion usually breaks down in the immediate vicinity of tree-level 
phase-space boundaries. Practically, this means that more work (read resumma- 
tions) is needed near the tree-level threshold {Q^ — AM^ ) in order to calculate 
accurately these processes. 

Concerning phenomenology, we find again that the off-shell annihilation is 
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dominant for very large dilepton energies. This new contribution therefore en- 
hances the thermal dilepton rates at moderate invariant masses and large ener- 
gies, and it would be interesting to include it in hydrodynamical evolution codes 
in order to see whether it leads to visible effects in a realistic heavy-ion collision 
scenario. 
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A Integration over p\ 



This is the simplest case since the complex conjugation in Eq. ( Pq) merely brings 
a factor 2. Indeed, the quantities u-f Af^g and (u-f w-t-M^j^)'^ — 4uu never vanish 
for positive u and v. In addition, we always have u + v + M^g > 0. In this case, 
the only nontrivial integral we need to perform is 



This integral is elementary. First rewrite (u + v + M^^)^ — Auv = (u — v + 
M^g-)^ + AvM^g. Then, introduce a new variable t defined by 



A.l Case > 




(65) 




(66) 



which brings the integral A to the form 




(67) 



at which point the integration is trivial: 



At this stage, it is sufficient to use the relation 




(68) 




(69) 



in order to obtain the first of Eqs. (^7|). 
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A.2 Case M^j < 

This case is more complicated that the previous one, mainly because A = (u + 
V + Mgfj)2 — Auv and u + v + M^g can both become negative. Let us first note 
that 1 / should in fact be interpreted asQ 

-^ = ^-.s.gn(po.„)^. (70) 



iTTsign{poro)Siu + M^ff) . (71) 



Similarly, we must write: 
1 



u + Af 2^ u + Af, 



2 

off 



Therefore, we have for the integral A the following expression 



A^2 du 



.(A)sign(. + . + A.2^) V Oi^ 



(72) 



One notes that in this equation, the second term corresponds to the cuts denoted 
(a) and (c) in figure |[ This is the reason why this term is always zero if 
Q2 < AMI. 

A.2.1 u + 4Af2fj>0 

Let us note first that this integral simplifies if w > —4M^g since then A and 
u+v+M^ff are both positive at all u. Using now the following change of variable 

and employing a strategy similar to the one used before, we find easily 



A„+4Af^„>o = , , tanh-^/ '^'^^^'ff , (74) 

which proves the first of Eqs. (|3^). 
A.2.2 V + 4Af < 

Things are more involved if + 4Af < since now both A and u + v + M^g 
can become negative. The second term in Eq. (|7^) is now different from zero 
and is given by: 

A. = -2.4=^=^ . (75) 

^/-viv + AM-g) 



^^The sign in front of the imaginary part is controlled by the sign of the infinitesimal 
imaginary part of M^jj. 
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In order to calculate the first term of Eq. 
of A as well as the sign oi u + v + M^. 
[0, v-M^g -2y^-vM'^f^] orue[v- 
2 ^/—vM^g, one can check that u + 
for u e [0,v- - 2V-wA/2ff] 
negative if u < —M'^^. For u ^ [v — 



we need to study the zeros 



We find that A is positive for u G 

oof 



P - 



For u > w - Af - 



M^g is always positive. On the contrary. 
'r V + M^ff is positive if u > 
2 



M^g and 



/f2 



-vM^g, +oo[, the appropriate 
change of variable is again the one given in Eq. ([73|) , which leads to the following 
contribution to A: 



+00 



du 



1 



" + ^^cff [iu + v + A.f2^)2 - 4ww]^/^ 



tan 



V - AM^g 



For we [0, w — Af2^ — 2^y—vM^g], one must define instead 



MO -cosh- 

By the same method, one can check that 



(76) 



(77) 



I 



du 



sign(M + V + M^g) 



^* + ^cff [{u + V + A'Qgf - Auv] 



1/2 



Stan ^ 

2 



(78) 



Adding up the three contributions of Eqs. (|75|), ( [76| ) and (fz^), we obtain the 
following result for A in the region where v + 4:M^g < 0: 



A, 



^~v{v + ur^g) 



itan 



-V - Ul^g 



(79) 



which proves the second of Eqs. (38). It is instructive to note that we have a 
cancellation between the first and second term of Eq. ( fl% ) of a term that would 
have been singular in the limit of zero momentum transfer {v — > 0). This is 
nothing but a manifestation of the KLN theorem for the thermal production 
of a massive particle, since these two terms correspond respectively to real and 
virtual cuts. Therefore, this is yet another example to fuel the controversy 
between |^,^ 65 and Q , which does not support the finding by [ |63| , |64p5[ | that 
the Kinoshita-Lee-Nauenberg |^,|l| theorem is invalid. Note that the present 
verification of the KLN theorem, for bremsstrahlung and off-shell annihilation 
processes, involves diagrams that in fact appear for the first time at three loops 
in the bare perturbative expansion. 
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